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Abstract 

We construct nonlinear coherent states for the Susskind-Glogower operators by the application of the 
displacement operator on the vacuum state. We also construct nonlinear coherent states as eigenfunctions of 
a Hamiltonian constructed with the Susskind-Glogower operators. We generalize the solution of the eigen- 
function problem to an arbitrary |m) initial condition. To analyze the obtained results, we plot the Husimi 
Q function, the photon number probability distribution and the Mandel Q-parameter. For both cases, we 
find that the constructed states exhibit interesting nonclassical features, such as amplitude squeezing and 
quantum interferences due to a self-splitting into two coherent-like states. Additionally, we show that non- 
linear coherent states may be modeled by propagating light in semi-infinite arrays of optical fibers. 

Keywords: Nonlinear coherent states; Susskind-Glogower operators; quantum squeezing; quantum inter- 
ference. 



1 Introduction 

Over the years, major effort has been directed towards the generation of nonclassical states of electromagnetic 
fields, in which certain observables exhibit less fluctuations than in a coherent state, whose noise is referred to 
as the standard quantum limit (SQL). Nonclassical states that have attracted the greatest interest include (a) 
macroscopic quantum superpositions of quasiclassical coherent states with different mean phases or amplitudes, 
also called " Schrodinger cats" [3 [2] , (b) squeezed states [5] , whose fluctuations in one quadrature or the ampli- 
tude are reduced beyond the SQL, and (c) the particularly important limit of extreme squeezing, i.e. Fock or 
number states [4]. 

In this frame, the coherent states for the electromagnetic field, introduced by Glauber [Fl IB] and Sudarshan 
[7], have been very important. Besides, their importance is increased by the fact that these states are relatively 
easy to produce in the laboratory and in their classical wave behavior. These states can be obtained by different 
mathematical definitions: 

a) As the right-hand eigenstates of the annihilation operator [B] 



with a a complex number. 

b) As the those states obtained by the application of the displacement operator D (a) = e""*"" ° on the vacuum 
state of the harmonic oscillator [B] 

\a)=Dia)\0). 



c) As those states which their time-evolving wave function shape does not change with time and whose centroid 
follows the motion of a classical point particle in a harmonic oscillator potential [8]. 
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Using the harmonic oscillator algebra, one can obtain the same coherent states from the three different 
definitions; however, for systems with complex dynamical properties, the harmonic oscillator model is not 
adequate, therefore new methods to generalize the idea of coherent states for systems like those have been 
proposed. 

Nieto and Simmons [9l[T0l|Tl] constructed coherent states for potentials whose energy spectra have unequally 
spaced energy levels, such as the Poschl- Teller potential, the harmonic oscillator with centripetal barrier and 
the Morse potential. 

Gazeau and Klauder [12j proposed a generalization for systems with one degree of freedom possessing discrete 
and continuous spectra. These states were constructed performing a parametrization of the coherent states by 
two real values: an amplitude J, and a phase 7, instead of a complex value a . 

Man'ko et aL[T3] introduced coherent states of an /-deformed algebra as eigenstates of the annihilation 
operator A = af (n) where n = a^a is the number operator and a, are the annihilation and creation boson 
operators of the harmonic oscillator algebra, respectively. A remarkable result is that these states present 
nonclassical properties such as squeezing and antibunching |14) . 

Recamier et al. |15) proposed coherent states using a deformed version of the displacement operator method 
generalized to the case of /-deformed oscillators, assuming that the number operator function appearing in the 
commutator between the deformed operators can be replaced by a number. The method yields a displacement 
operator which is approximately unitary and displaces the deformed annihilation A and creation A^ operators 
in the usual way. 

In this work, we construct nonlinear coherent states by the application of the displacement operator (for 
the Susskind-Glogower operators [16]) on the vacuum state and time-dependent |m) displaced number states 
as eigenfunctions of a Hamiltonian representing the fundamental physical coupling to the radiation field via 
the Susskind-Glogower operators. We show the Husimi Q function [17] for the resulting states, as well as 
their photon distribution and Mandel Q-parameter [18] , in order to determine the nonclassical features of the 
constructed states. Additionally, we show that nonlinear coherent states may be modeled by propagating light 
in semi-infinite arrays of optical fibers. 

Nonlinear coherent states may be constructed using the Susskind-Glogower (SG) operators, for instance, 
defining a displacement operator for them acting on the vacuum state. The Susskind and Glogower proposed 
operators are 



00 



V^y^n){n + l\^^==a, (1) 
V"- + 1 



n=0 



= y\n+l){n\^a)^==, (2) 



=0 



satisfying the conditions 



V\n) = \n-l), (3) 
V^\n) = \n + l). (4) 
Additionally, we would like to make explicit the result 

V\0)=Q, (5) 

that comes naturally from ([1]). 

The SG operators possess a non-commuting and non-unitary nature, that resides in the expressions 

VV^ = 1, (6) 

y^y=l-|0)(0|. (7) 

From the above expressions we can see that the non-commuting and non- unitary nature of SG operators is 
only apparent for states of the radiation field that have a significant overlap with the vacuum 

= (^|0)(0|V'). (8) 

Therefore, for states where the vacuum contribution is negligible, we can consider them as unitary and commu- 
tative and we can perform the following approximation 

v-^ ~ VK (9) 
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The properties of the SG operators play an important role in the development of the present work. For 
instance, if we analyze ([1]) and we find that SG operators have the same form of the ones we need to 
construct nonlinear coherent states, i.e., V ^ f {h + 1) a and V'' = a^/ {n + 1) . 

Following Recamier et a/.[TS] we define SG coherent states (from this point, we will refer to these states as 
SG coherent states, keeping in mind that, in fact, they present a nonlinear behavior) by 

\a),^ = e<^'-^)\0), xeR. (10) 

As we have seen, commutation relations for the SG operators are not simple, so we cannot factorize the 
displacement operator in a simple way, which is why we propose two methods to solve the displacement operator. 
First, we use the approximation V^^ ~ to solve, approximately, the displacement operator. This solution 
helps us to understand how the exact solution for the displacement operator should be. Second, we solve the 
displacement operator in an exact way by developing it in a Taylor series, that allows us to introduce the exact 
solution for nonlinear coherent states, constructed with the SG operators. Finally, we analyze the constructed 
states via the Q function |17| . the photon number distribution and the Mandel Q-parameter [18| in order to 
show their nonclassical features such as amplitude squeezing and quantum interferences. 



2 Approximated displacement operator 

A first approach to factorize the displacement operator in the product of exponentials, is to consider the 

approximation V^^ ~ V'' . 

Let us write the displacement operator as 

^,,^eK^-M-^ (n) 

We find that the right-hand side of this equation corresponds to the generating function of Bessel functions, 
that implies that 

oo 

DsG- V^"Jn{'2x), (12) 

n— — oo 

where Jn is the Bessel function of the first kind and order n. 
Applying the displacement operator on the vacuum state, we have 

oo 

i?SG|0)~co J2 V^"Jn{2x)\0}, (13) 

n— — OO 

and using that V'^ \n) = |n + 1) and that V"-^ V\ we obtain 

oo 

|a)sG=^SG|0)^co^ J„(2a;)|n). (14) 
From the normalization requirement we determine cq. 



^ J„ (2x) (n| J2 Jrn (2x) |m) 

n— m— 

oc 

= clY,Jl{2x)^l. (15) 



n=0 



Using the result [19] 



we obtain 



1 = j2(2a;) + 2^ J2(22,)^ (16) 
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and substituting in (HH), we finally get 



As we mentioned before, solution (|18p helps us to foresee the exact solution for these states. We can expect 
that the exact solution for SG coherent states corresponds to a linear combination of number states where the 
coefficients are, except for some terms, Bessel functions of the first kind and order n. 



3 Exact solution for the displacement operator 

In order to "disentangle" the displacement operator in an exact way, we can develop the exponential ()10|) in a 
Taylor series and then to evaluate the terms ^V^''' — . For instance, for = 7 we have 



7 /7 



1) (0|-|0>(1|) 

(19) 

;j(|3) (0|-|2)(1| + |1)(2|-|0) (3|) 
+ Q (|5) (0| - |4) (1| + |3) (2| - |2) (3| + |1) (4| - |0) (5|), 

where : : means to arrange terms in such a way that the powers of the operator V are always at the left of the 
powers of the operator . 

From the definition of the SG coherent states (fTO|). we can write 

la; 



SG — 



k=0 n=0 



(20) 



kl 

k=0 n=0 



where the square brackets in the sum stand for the floor function, which maps a real number to the largest 
previous integer. 

We can rewrite the above equation as 



- -i ^ r'^ ^ /k\ - 1 (21) 



fc=0 ' n=0 ^ 



where we have taken the second sum to oo as we would add only zeros. We now exchange the order of the sums 

I")SG = 

= e-'^^e^^' |0) +V'yy -tll^V^-V^" |0) . ^^^^ 
' ' ^ ^ (k- n)\n\ ' ' 

n=0 k=n ^ ' 

By setting m = k — n and using that 

T>2«l/t'"+" |o) ^ f t" |o) , (23) 



4 



we get 



/SG — 

oo oo 

-xV xV^ 



^-^ ^-^ m'.nl 

n— m— 

e--e-^ 10) + E 10) 



n! — ' m! 

n— m— 



-xv^xv^ |0) + T/^e-^^e"^' |0) , 



to finally obtain 



\a)sG={^ + V')e-^^e^^'\0). (24) 
Applying the exponential terms on the vacuum state, we have 



^)sG = (l + V') E (2a;) \n) 

oo oo 

= ^ J„ (2x) \n) + E (22;) h - 2) , 



(25) 



n=0 n=2 

making to = n — 2 in the second sum and performing the index change n = k ~ 1, we obtain 

oo 

l")sG = E t'^'^-i (2^) + -^fc+i (2a;)] Ifc - 1) ■ (26) 

k=l 

Using the recurrence relation of the Bessel functions [12] 

xJn-i (x) + xJn+i {x) = 2nJ„ (x) , (27) 
and changing again the summation index, we finally write 

^ oo 

l")sG = -E("+l)-^«+i(22;)l")- (28) 

Equation (j28p is an important result because it constitutes a new expression for nonlinear coherent states. It 
remains to analyze the behavior of the constructed states in order to determine the nonclassical features that 
nonlinear coherent states may exhibit. 

Before we proceed with the analysis, and since we will need this result later, we make clear that, as we can 
verify from ([28]) . for a; = 0, we have 

l«(^ = 0))sG = |0). (29) 



4 SG coherent states analysis 

There are different ways to find out if the state we are constructing resembles one that we already know. Here, 
we will use three different methods: the Husimi Q function [17| . the photon number distribution and the Mandel 
Q-parameter [T5] . 



4.1 The Q function 

The Q function, introduced by Husimi [17], corresponds to a quasiprobability function that helps us to determine 
the behavior of a quantum state in phase space. The Q function is defined as the coherent state expectation 
value of the density operator and is given by 

Q^- {a\p\a) . (30) 
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t 




Fi gurc ll The Q function for (a) coherent state |ci) with a = 2 and (b) number state \n) with n = 5. 



We plot the Q function for a coherent state |a) and for a number state \n) in Figure [TJ If we substitute 
p — \ip) for the SG coherent states in ([5D|) . we obtain 



QsG {a,t) = 



(31) 



Figure m shows the SG coherent states Q function for different values of the parameter x. 

We observe that the initial coherent state squeezes. Later, we will find the value of x for which we obtain 
the maximum squeezing of the coherent state. We can also see that, as the parameter x gets bigger, the state 
splits into two coherent-like states, this produces quantum interferences, as we will see later. 



4.2 Photon number distribution 

When one studies a quantum state, it is important to know about its photon statistics. The photon number 
probability distribution P {n) ~ {n\a){a\n) ~ |(n|Q;)|^ is useful to determine amplitude squeezing. We should 
refer to amplitude squeezed light as light for which the photon number distribution is usually narrower than 
the one of a coherent state of the same amplitude. The photon number distribution is also useful to analyze 



if there exist effects due to quantum interferences. Using P (n) 
photon number distribution as 

P{n)= - (n + 1) Jn+i {2x) 



|(n|Q:)| , we write the SG coherent states 



(32) 



Figure |3] shows the SG coherent states photon number distribution for different values of the amplitude param- 
eter X. 

Figure [3] helps us to understand the effect of quantum interferences; for instance, consider Figure [Sjc) . We 
see that it is not a uniform distribution of photons, the distribution has "holes" ; these holes are the consequence 
of the interference between the two states arising from the splitting of the initial one. 

Comparing Figure |3l[a) with the one obtained for a coherent state wave function that is moving through 
the harmonic oscillator potential between the classical turning points, we can see that the photon number 
distribution for the SG coherent states is narrower than the one for a coherent state of the same amplitude. 
This suggests that we are, indeed, obtaining an amplitude squeezed state. It is interesting to know when the 
state is maximally squeezed, but we need another tool to obtain the value of the parameter x for which this 
occurs. 
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Figure 2: Exact SG coherent states Q function for (a) x = 1; (h) x = 5; (c) a; = 10 and (d) x = 20. 



4.3 Mandel Q-parameter 

There has been an extensive argument about which is the better way to determine the quantumness of a given 
state; nevertheless, there exists a very useful tool for determining the nature of the states we have constructed. 
This tool is called the Mandel Q-parameter [18] . We will use it not only because it represents a good parameter 
to define the quantumness of SG coherent states but also, because it will allow us to find the domain of x for 
which they exhibit a nonclassical behavior; moreover, it will help us to find the value of x for which the state 
is maximally squeezed. 
The Mandel Q-parameter is defined by 

Q = ^ - 1' (33) 



where 



If Q 



>0, 
-0, 
<0, 



-1, 



super Poissonian distribution 
Poissonian distribution (coherent state) 
sub-Poissonian 
number state. 



(34) 



For the SG coherent states we have that 



Lfc=0 



-Y^k^Jl(2x) 



k=0 



(35) 
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Figure 3: SG coherent states photon number probabihty distributions for (a) x = 1; (b) x = 5; (c) x = 10 and (d) x = 20. 



and 



£ fcV,2 (2x) - 2 £ k'ji {2x) + £ (2a 



Lfc=0 



fe=0 



fc=0 



(36) 



The even sums, with respect to the power of fc, in ((35|) and ([36l) . can be evaluated. 
In the Appendix, we show explicitly that 



k=0 



and 



^ k^Ji i2x) = 3x^ + 



(37) 



(38) 



The sum X^fc^o ^^"^fc ('^•^) more complicated, and can be evaluated using the technics developed by DattoUi 
et al Eni, to 



k^Jl {2x) = a;^{(6a;2 + 1) Jo (2a;) + (G-x^ - 1) Ji (2a;) 



k=l 



(39) 



2x' 



2xJo(2x)Ji(2x) + — [Jo(2a;)J2(2a;) + Ji(2a;) J3(2a;)]}. 
o 



Substituting the values of the sums into equation (|33)). we obtain the plot shown in Figure ID From Figure |4] we 
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Figure 4: Mandel Q-parameter for the SG coherent states. 



can see that, depending on the parameter a;, the photon distribution of the constructed states is sub-poissonian, 
Q < 0, meaning that amphtude squeezing states may be find for a value of x within the domain < a; < 13.48. 
Also, we find that the most squeezed state may be obtained at a; = 2.32, with Q — —0.64. 
At this point, all the comments on the results have been made considering x as a parameter; however, in what 
follows, we will see that x corresponds to an expression that depends explicitly on what we should refer to as 
interaction time. 

5 Eigenfunctions of the SG Hamiltonian 

Previously we managed to construct nonlinear coherent states applying the displacement operator on the vac- 
uum state. However, even when the obtained results are very interesting, we cannot avoid to wonder how these 
states could be physically interpreted. Physical interpretations are given by operators representing observables, 
i.e., quantities that can be measured in the laboratory. The most important observable is the Hamiltonian, this 
operator helps us to find the energy distribution of an state via its eigenvalues. 

Here, we construct SG coherent states as eigenfunctions of a Hamiltonian that we propose and which 
represents the fundamental coupling to the radiation field via the SG operators. In section 5.1 we construct 
time-dependent SG coherent states with the vacuum state |0) as initial condition, i.e., we construct states that 
satisfy |a (x = 0))gQ = |0). In section 5.2 we generalize the eigenfunction problem for an arbitrary \m) initial 
condition, we also show that previous results correspond to the particular case m = and in section 5.3 we make 
use of the three methods presented in the previous chapter in order to analyze the properties and time-evolution 
of the constructed states. 

5.1 Solution for |0) as initial condition 

As we mentioned before, it is possible to construct SG coherent states as eigenfunctions of the interaction 
Hamiltonian 



where is the coupling coefficient. Hamiltonians like this may be produced in ion-traps [STJ [551 [53] . 

The Hamiltonian proposed in (j^O)) corresponds to a variation of the one used in [53] to model physical couplings 

to the radiation field. Here, physical couplings take place via the SG operators. 




(40) 
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We write the eigenfunctions of the Hamihonian (PD|) in the interaction picture as 

oo 

= (41) 

Tl = 

and we have 

oo oo 

^|^)=77^C„|7i-l)+r;^C„|n+l). (42) 

n— 1 n—0 

Now, changing the summation indexes, we obtain 



(43) 



(44) 



F IV') = ^ C„+i |n) +tjY1 ^n-i \n) 

n—0 n— 1 

oo 

= r,Ci |0) + ?7 ^ {Cn+i + C„-i) \n) . 

n=l 

Then 

OO 

7lCi\Q)+Y,V{Cn+l+Cn-l)\n) = 
n^l 

oo 

= ECo\0} + Y,ECn\n}. 

Comparing coefficients with same number states of the sum, we have 

Ci = -Co, (45) 

ry(C„+i + C„^i) = SC„. (46) 
These are the recurrence relation of the Chebyshev polynomials of the second kind (|19)), and we can write 

oo 

|V'(i;C))-E^~'''*^"(^)l")' (47) 

n=0 

where 

E = 2T]t (48) 

However (|47| does not satisfies the initial condition |-0 (i = 0)) — |0). Moreover, the solution (j47|) has the 
parameter ^ and, as we see from (jlOp . we should not have another parameter, except for the time t. 

A way to construct a solution, as the one we previously obtained in (|28p. is by looking at the exponential 
term in (j47p and noticing that it has the form of the Fourier transform kernel; so, we need to propose a ^- 
dependent function and integrate it over all ^, in order to obtain a solution where time is the only variable. We 
have then, 

/-oo 

IV [t)) '^^P (^) ^« (^) I") • (49) 

J-OO 

We see from the above equation that |?/; [t)) corresponds to a sum of Fourier transforms of Chebyshev polynomials 
with respect to a weight function P (^). This kind of Fourier transforms may be solved by using the following 
result [21 



Jn (w) 



= J^- HT Un-i (0 yr^rect f I 

V TT 71 \Z 



(50) 
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where is the Fourier transform and 



'^'^ct [ 2 j jo otherwise. ^'^^'^ 



From (|50l) we write 



J„ (lo) 



1 



(52) 



Using definition (1511) , and making k ~ n ~ I, we obtain 

Jk+i (w) 



1 ■ !■ .xfc+1 
1 



[/, (0 e^'^«de (53) 



TT fc + 1 

With equation ([53)) it is possible to solve the integral in (|49|. Writing the weight function as 



^(O-§yT^rect(^0 



(54) 



substituting it into equation (|49l) and using now ([53)) . we get 

2 



9 _ ^ 

1^ = E + 1) ■^"+1 (-2'?*) I") • (55) 

' n=0 



Considering the odd parity of the Bessel functions, we finally obtain 



1^ = - 5] (n + 1) J„+i (277t) |n) . (56) 
We see that \ip {t)) in equation ([56]) depends only on t and considering that 

\ i , „ = 1, 

we can verify that 

\^{t = 0)) = \0). (58) 
The equation ([56]) corresponds to the expression for SG coherent states that we obtained previously ([28]) . 

The solution presented in this section allows us to notice that, while in the previous chapter x was only a 
parameter, now it represents something physical. It may be related to an interaction time, for example, in the 
motion of a trapped atom [53]. We have managed to construct the same expression for the SG coherent states 
as the one obtained by the application of the displacement operator on the vacuum state; however, we will 
see that the formalism presented in this section may be used to generalize the solution for an arbitrary initial 
condition |m). 

5.2 Solution for |m) as initial condition 

At this point, we have managed to construct SG coherent states, first as those obtained by the application of 
the displacement operator on the vacuum state and later, as eigenfunctions of the Hamiltonian ([10]) : however, 
they are a particular case of a more general expression. 
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Using the recurrence relation for the Chebyshev polynomials and the result (|55|) . it is possible to generalize SG 
coherent states to an arbitrary |m) initial condition, where to = 0, 1, 2, .... 
From equation (jlH]), we have that 

IV' (t = 0)) = ^ / d^P iO C/„ iO \n) , (59) 
„=o"'-°° 

and considering the function 



Prn iO = lV^~~eU,n {0 TCCt (^^ 



= E /°° -Vl^Uo (0 f/« (0 e-^^"«* rect f |) I 
= E / - VT^C/o (?) C/n (C) e-'^^^^dC |n) 

— J 1 TT 



n=0 



For TO = 1, we have 



= E/ PiiOUniOe-'^'^^'d^ln) 

= E -yr^C/i (0 Un (0 e-*"'^* rect (|) \n) 

= E /' -yT^(20t^n(C)e-*'"«*dC|n). 



(60) 



and the well known Chebyshev polynomials of the second kind orthonormal condition 

- c/n (0 (0 yr^rfe = '5„„, 

we obtain 

oo 

\'4> {t = 0)) = E I") • (61) 

To obtain the solution for the m-state, let us consider the particular cases m = and m = 1. 
For TO = 0, 

IV'(i))™.o = 

oo /.oo 



(62) 



n=0 

Using (j52p and making n — > n + 1 , we have 

1^ (0)™=o = E ^" (n + 1) :^!^±liM |,) . (63) 

n=0 

Using the recurrence relation 

xJn-i (22:) + xJn+i (2x) = nj„ (2a;) , (64) 
and changing the summation index, we finally obtain 

oo 

IV' (0)™=o = E t*"-^" (2r;i) + J„+2 (2ryt)] |n) . (65) 



(66) 
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Using the recurrence relation for the Chebyshev polynomials of the second kind, we write 

CO p\ 



^ Pi cy 



Considering (|52)) for the first integral, and making n = n + 2 for the second one, we have 

= Y.-, [*"^' + 2) ^"+2 (2r?i) + z"-^nJ„ (2ryi)] |n) . 
We now construct a different expression for the above equation ([55]). Let us write it in the 

lV'(t))„.i = 



E 

n=0 



77< 



making /c = n + 2 



E 

n=0 



(k - 2) Jfe_2 (2r?i) + i 



fe-i 



77t 



fcJfe (2??i) 



using the recurrence relation (|64p , we write 

OO 

i^w)™=i = E 



(fc - 2) Jfc_2 (2r/t) 
[Jfc-i (277i) + Jfe+i (27?^)] 



,-fe-i j_ 



g| -^'=-^[i(fc-2)J.-2( 

E ^ 



= E 

Making n = A; — 2, we obtain 

OO 

E 



(2r7t) - Jfc_i (2r7<) 
(277t) 

(fc - 2) Jfc_2 {2vt) - Jfc-2+i (2?7i) 



+1^-2+1 Jfc_2+3(277t) 



rjt 



nJn {2r]t) - Jn+i (2?7<) 



Using again the recurrence relation (1641) , we finally write 

OO 

H {t))m=l = E (2^<) + *"+'^«+3 (2r?i)] 

Following the same procedure, it is easy to prove that for m = 2, 

OO 



n=0 
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Rewriting results (US]), ^ and dTJ), 

oo 

= E [^"""-^"-o (2r?i) + J„+o+2 (277^)] \n) , 

n=0 
I^W)™=1 = 

OO 



(75) 



n=0 



m=2 

oo 



n=0 



It is easy to see that the solution for the m-initial condition is 

= [i"~"'Jn-m {2vt) + J„+™+2 {2f^t)] \n) . ^^^^ 

We have constructed a new expression for nonlinear coherent states and that we call SG coherent states. We 
have managed to construct an expression that allows us to study the time evolution of SG coherent states for an 
arbitrary \m) initial condition. We also found the physical interpretation of the parameter x (used in previous 
sections) as a normalized interaction time with respect to the coupling strength, i.e., x = rjt. We analyze now 
the nonclassical features of the constructed states, so we have to make use of the methods previously mentioned, 
these are the Q function, the photon number distribution and the Mandel Q-parameter. 



5.3 Time-dependent SG coherent states analysis 

To perform a complete description of the constructed states (j76|) , we have to verify if they present the nonclassical 
features that nonlinear coherent states may exhibit. In order to study these nonclassical features, we propose 
to use three methods. First, we analyze their behavior in phase space via the Q function; then, because we 
want to analyze amplitude squeezing and quantum interferences, we show the photon number distribution of 
the constructed states, and finally, as we want to know when the constructed states are maximally squeezed, 
we show the Mandel Q-parameter. 



5.3.1 Q function 

Considering the definition Q — ^ (a|p|Q;) and writing p in terms of the constructed states (|76p . we write the Q 
function for the time-dependent SG coherent states as 

QsG (a) = 

°° *" ^ (77) 

V ^ J„-™ {2vt) + J„+,„+2 (277t)] . 

n=0 

Figure [5] shows the time-evolved Q function of SG coherent states for different initial conditions. The figure 
is structured as follows: The time evolution is shown in each row, for example, (a) represents the time evolution 
of SG coherent states for |0) as initial condition and the index (i) represents the normalized interaction time rjt. 
We have chosen (ii) to be iji = 2.32, because, as we showed previously, it is the time when the state with initial 
condition |0) is maximally squeezed. The nonclassical features of the constructed states are summarized in Fig- 
ure[S] SG coherent states present a strong amplitude squeezing (Figure [Sfa,ii)), splitting into two coherent-like 
states (Figure[Sl[b,ih)) and, as a consequence of the splitting, pronounced quantum interferences (Figure [SJc,iii)). 
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Figure 5: SG coherent states Q function for (i) rjt = 1; (ii) ryi = 2.32; (iii) rjt = 5 and (iv) rjt = 20, with initial conditions (a) 
(b) |1>; (c) |5> and (d) |10>. 



An interesting result is that no matter what initial condition we choose, SG coherent states eventually split 
into two coherent-like states. This is very interesting result because such a distribution corresponds to states 
called Schrodinger's cat states, which are very useful in quantum information processing [25] . 

5.3.2 Photon number distribution 

To complement the description of the nonclassical features that we observed from the Q function, we show in 
Figure [5] the photon number distribution of the SG coherent states considering the same conditions of Figure 

m 

Time-dependent SG coherent states photon number distribution for the to- initial condition is given by 

Pm {n,t) = |z"-'"J„_™ (277t) + i"+'"J„+„+2 (2r7i)f , (78) 



5.3.3 Mandel Q-parameter 

As we want to know the time domain for which the constructed states exhibit amplitude squeezing, and moreover, 
we want to know when the states are maximally squeezed, we obtain the Mandel Q-paxametev for the time- 
dependent SG coherent states. 
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Fi gure 6; SG coherent states photon number probability distributions for (i) rjt = 1; (ii) rft = 2.32; (iii) r}t = 5 and (iv) rjt = 20, 
with initial conditions (a) |0); (b) |1); (c) |5) and (d) |10). 



We have that 

oo 

(n) = ^ n [Jn^ra {2Tjt) + (-1)'" J„+™+2 {2Tjt)f , (79) 

n=0 

and 

oo 

= E ["^"-™ (2^70 + (-1)" Jn+m+2 {2vt)f . (80) 

Substituting in the definition of tlie Mandel Q-parameter, equation ([33]), we obtain the plot shown in Figure [Tj 
From Figure[7]we see that (a) shows the Q-parameter for the SG coherent states where the maximum squeezing 
happens when rjt — 2.32. From the others, we cannot observe squeezing, because the initial condition is a full 
squeezed state, a number state; however, we obtain that an initial number state eventually transforms into 
two coherent-like states that we may identify as Schrodinger's cat states and, due to the splitting, we observe 
quantum interferences. 

The generalization (j76l) does not give us different effects from the ones that we may obtain from (1561) : however, 
as we will see in the next section, equation ()76|) helps us to show that nonlinear coherent states may be modeled 
by propagating light in semi-infinite arrays of optical fibers. 
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Fi gure 7'. SG states Mandcl Q-parameter with initial conditions (a) |0); (b) |1); (c) |5) and (d) |10). 



6 Classical quantum analogies 

The modeling of quantum mechanical systems with classical optics is a topic that has attracted interest recently. 
Along these lines Man'ko et al. [26| have proposed to realize quantum computation by quantum like systems, 
Chavez-Cerda et al [27] have shown how quantum-like entanglement may be realized in classical optics, and 
Crasser et a/. [28] have pointed out the similarities between quantum mechanics and Fresnel optics in phase 
space. Following these cross-applications, here we show that nonlinear coherent states may be modeled by 
propagating light in semi-infinite arrays of optical fibers. 

Makris et aZ.|29| have shown that for a semi-infinite array of optical fibers, the normalized modal amplitude in 
the nth optical fiber (after the mth has been initially excited) is written as 

an {Z) = Ao [i"-" J„_™ (2Z) -t- r+™ J„+„+2 (2Z)] , (81) 

where Z — cz is the normalized propagation distance with respect to the coupling coefficient c. 
We see that, for Aq ^ 1 the normalized intensity distribution is 

/„ (Z) - |z"-"J„_„ (2Z) + z"+"J„+™+2 (2Z)|' . (82) 

Rewriting (j78p and using the normalized interaction time with respect to the coupling coefficient ry, i.e., x — r]t, 
we have 

Pm{n,x) = |z"-™J„_™(2x)+z"+™J„+„+2(2a;)|'. (83) 

As equations (|82|) and (|83|) are the same, we conclude that the photon number distribution for the SG coherent 
states may be modeled by the intensity distribution of propagating light in semi-infinite arrays of optical fibers. 
We have found a new relation between quantum mechanical systems and classical optics. 

7 Conclusions 

We managed to construct a new expression for nonlinear coherent states, that we called SG coherent states, by 
developing the displacement operator in a Taylor series. The importance of nonlinear coherent states resides in 
the nonclassical features that they may exhibit and, in order to analyze the nonclassical behavior of the con- 
structed states, we presented three methods: the Q function, the photon number distribution and the Mandel 
Q-parameter. 

The expression for the SG coherent states depends on a parameter x and, considering the values of a;, we 
found that 
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• For a certain doraain of x, SG coherent states exhibit amphtude squeezing; i.e., their photon number 
distribution is narrower than the one of a coherent state of the same amphtude. Using the Mandel 
Q-parameter we determined that maximum squeezing occurs when x = 2.32. 

• For larger values of x, SG coherent states split into two coherent-like states, this superposition of coherent 
states is called Schrodinger's cat states and are very useful in quantum information processing [25] . We 
also found that the splitting of the states gives rise to quantum interference effects. 

We also constructed SG coherent states as eigenfunctions of a Hamiltonian for the SG operators represent- 
ing physical couplings to the radiation field. We managed to obtain the solution as the one obtained before, 
except for a phase term that rotates by 7r/2 the Q function, but it does not change the results that we have 
already found. The obtained expression for SG coherent states helped us to physically interpret the parameter 
a; as a normalized interaction time with respect to the coupling coefficient 77. 

The formalism presented in this part allowed us to construct a solution for an arbitrary initial condition \m). 
Equation (j76|l is the most important result of this work, because it represents a new expression for nonlinear 
coherent states; moreover, it represents a general expression for an arbitrary |m) initial condition. 

Analyzing expression (j76p . we found that a multiple splitting of the states occurs; this gives rise to interesting 
structures of the phase-space distribution as the one shown in Figure [5ljc,iii). Except for that result, the 
generalized solution presented the same effects as the particular case m = 0; however, it turns out to be 
very useful to establish an interesting relation between quantum and classical optics. We showed that SG 
coherent states photon number distribution may be modeled by the intensity distribution of propagating light 
in semi-infinite arrays of optical fibers. With this finding, we have presented a new analogy between quantum 
mechanical systems and classical optics. 

Appendix 

Sums of the Bessel functions of the first kind of integer order 

We will derive in this appendix the solution of some sums of Bessel functions of the first kind of integer order 
that appear in several applications, and in particular, that appear in this contribution in the section where we 
calculate the Mandel Q— parameter (section 4.3). We will demonstrate that 



where j/ is a positive integer, g (y) — ixsiny and Bn (xi,X2, ■■■,Xn) is the complete Bell polynomial [501 \ST[ [5^ 
given by the following determinant: 



00 



k=l 




(84) 



— TT 
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Bn{xi,X2, ■■■,Xn) 



= det 





r,')x2 




("-1)X4 




-1 


X\ 











-1 


X\ 


r,')x. 
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-1 Xl 



(85) 



To demonstrate (|84| , we will need the well known Jacobi- Anger expansions for the Bessel functions of the 
first kind ([33], page 933, [34], page 361; [19], page 70), 

oo 

e'^'°^y ^ i"J„(x)e"'' (86) 



and 

Using expression ([86]) . we can easily write 



^ Jn{x)e'"''. (87) 



n— — oo 



To calculate the n-derivative in the left side of equation above, we use the Faa di Bruno's formula ([33], page 
22) for the n-derivative of the composition 

, s (89) 

= E /'"^ (5 i^)) ■ Br^^k (<?' (x) , g" {x) , (x)) , 

fc=0 

where Bn^k ixi,X2, ...,a;„_fc+i) is a Bell polynomial [501 ISIl 131] ) given by 

Bn,k {xi,X2, a;,i-fe+i) = 



ji!i2!...jn-fc+i! (90) 
Ml!) (21) 



'Xiyi/XaV^ f Xn-k+1 

(n-fc + 1)! 



the sum extending over all sequences ji, ja, js, jn-fc+i of non- negative integers such that ji+j2 + ---+jn-k+i 
k and ji + 2j2 + 3j3 + ... + (n - + 1) j„_fe+i = n. 
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Using 



^ ^ix sin y 

(91) 

= e"^"^y ^ (.g' (x) , g" (x) , g^"-'^+^'> (x)) . 

k=Q 

We multiply now for the complex conjugate of (j87p and obtain 

oo 

fe,/=-oo (92) 
= S„ (g' (y),g" (2/),...,5(") (y)) . 

Integrating both sides of the above equation from — tt to tt, and using that e*('^~'^^(i?/ = (5fei, we arrive to 
the formula we wanted 



TT 

= ^ / (ff' (y) , 5" (y) , -,9^"'^ (y)) dy. 



(93) 



In particular, as the complete Bell polynomials for n = 2 and n = 4, are 

B2{X1,X2) ^xl+X2 (94) 

and 

i34(a;i, a;2, X3, 0:4) — x^ + Qx\x2 + 4x1X3 + "ix^ + X4, (95) 

it is very easy to show that, 

00 _j 

Y^k'jUx)^^ (96) 



4 

fc=i 



and 



00 „ 

^fcV,^(x) = A^4 + ix^. (97) 
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